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Abstract. The state of the art in the study of n" , K~ and E~ atoms, along with the in- 
medium nuclear interactions deduced for these hadrons, is reviewed. A special emphasis is placed 
on recent developments in A'-nuclear physics, where a strongly attractive density dependent 
if _ -nuclear potential of order 150-200 MeV in nuclear matter emerges by fitting if~-atom 
data. This has interesting repercussions on K quasibound nuclear states, on the composition of 
strange hadronic matter and on K condensation in self bound hadronic systems. 



1. Introduction 



Hadronic atoms have played an important role in elucidating in-medium properties of hadron- 
nucleon (hN) interactions near threshold [1, 2]. Hadronic atom data consist primarily of strong 
interaction level shifts, widths and yields, derived from /i _ -atom X-ray transitions. These data 
are analyzed in terms of optical potentials Vj^ pt = tf l jq(p)p which are functionals of the nuclear 
density p(r), capable of handling large data sets across the periodic table in order to identify 
characteristic entities that may provide a link between experiment and microscopic approaches. 
Here, tf l j^(p) is a density dependent in- medium hN t matrix at threshold, satisfying the low- 
density limit thjsf(p) — > tfrx, with the free-space t matrix i^jy, as p — > 0. A schematic summary 
of lessons gained by analyzing the available hadronic atom data in terms of optical potentials 
is given in Table 1, where the number of data points included in these analyses is shown in 
parentheses in the last column. A more exhaustive discussion of 7r~, S _ and K~ atoms follows 
below. 



Table 1. Hadronic atom scenarios which except for p atoms are discussed in the present review. 
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2. Partial restoration of chiral symmetry from pionic atoms 




Density independent optical potential global fits to pionic atom data across the periodic table 
reveal an anomalous s-wave repulsion [1, 2], a major component of which is due to a too repulsive 
isovector ttN amplitude b\ with respect to the free-space value bi ee . This is demonstrated on 
the l.h.s. of Fig. 1 which shows results of global fits to pionic atom data as a function of a 
parameter 7, related to the difference between neutron and proton rms radii: 

N-Z . 

r n -r p = 7 — h 6. (1) 

Applying a finite-range folding with rms radius 0.9 fm to the p-wave part of V£L t , the resulting 
X 2 minimum for a 'skin' neutron distribution is obtained at 7 = 1.1 ± 0.1 fm. 1 The r.h.s. of 
the figure shows results of global fits with empirical energy dependence imposed on the s-wave 
amplitudes 60 and b\ [4], and more importantly with a DD renormalization of b\\ 

h(p)= , h \ p , (2) 
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where f n = 92.4 MeV is the pion weak decay constant and a ~ 50 MeV is the ttN a term. Eq. (2) 
was derived by Weise [5] considering an in-medium extension of the Tomozawa- Weinberg (TW) 
LO chiral limit for b\ in terms of f n [6] which is then related to the quark condensate < qq >: 

MP) = -FWT. 4^ = 1^.1-^. (3) 

8vr/2( p ) /2 <qq> mlfl 

The figure makes it evident that the magnitude of b\ on the r.h.s., following the DD 
renormalization, is systematically smaller than that on the l.h.s., and at the \ 2 minimum it 
agrees perfectly with b { { ee . 



1 For a recent discussion of the role of neutron distributions in hadronic atoms, see Ref. [3]. 



Table 2. Values of b\ derived from density- independent fits to pionic atom data sets listed in 
Ref. [9]. For comparison, 6^ rcc = -0.0864 ± 0.0010 m" 1 [10]. 'Deep' refers to Is vr~ 'deeply 
bound' states in 205 Pb [11] and 115 > 119 > 123 Sn [8] observed in (d, 3 He) experiments at GSI. 
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A similar conclusion was reached in Refs. [7, 8] from measurements of Is 'deeply bound' 
pionic atoms of Sn isotopes. The advantage of using Is 'deeply bound' levels is that the p-wave 
irN interaction plays there a secondary role, but this merit is more than compensated by the 
considearably increased errors associated with smaller data sets, as shown in Table 2. The 
uncertainty listed in the fourth column makes it clear that the 'deeply bound' atoms alone do 
not give conclusive evidence for the need to renormalize b\. In fact, Suzuki et al. [8] considered 
'deeply bound' Is levels in three Sn isotopes together with 'normal' Is levels in 16 0, 20 Ne and 
28 Si, 12 data points in total yielding b\ = — 0.1149±0.0074 m~ x . However, this small uncertainty 
excludes the uncertainty from the p-wave 7riV potential which was held fixed in their analysis. 
A more realistic uncertainty for this type of deduction is given in column 3 of the table. 
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Figure 2. Low energy tt^ nucleus elastic scattering angular distributions [12] reproduced with 
bi{p) ansatz, Eq. (2). 



The renormalization of b\ derived from pionic atoms is consistent with that shown recently 
to be also required in low energy ir- nucleus scattering, as demonstrated in Fig. 2 for 21.5 MeV 
7r ± scattered off isotopes of Si, Ca, Ni and Zr at PSI [12]. 

3. £ Nuclear repulsion from £ atoms 

A vast body of (K~,tt^) spectra indicate a repulsive and moderately absorptive £ nuclear 
potential V^, with a substantial isospin dependence [13, 14]. These data, including recent 
(tt~, K + ) spectra [15] and related DWIA analyses [16], provide credible evidence that £ hyperons 
generally do not bind in nuclei. A repulsive component of a DD £ nuclear potential was already 
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Figure 3. from fits to S atoms [2]. The transition from attraction to repulsion occurs well 
outside of R c , the half-density radius of the charge distribution in Ca (l.h.s.) and Pb (r.h.s.). 



deduced in the mid 1990s from X~ atom data [17, 18], as shown in Fig. 3. In fact, is 
attractive at low densities outside the nucleus, as enforced by the observed 'attractive' X - 
atomic level shifts, changing into repulsion on approach of the nuclear radius. The precise 
magnitude and shape of within the nucleus, however, are model dependent as demonstrated 
by the difference between potentials DD and F (defined in the Appendix). This repulsion 
bears interesting consequences for the balance of strangeness in the inner crust of neutron stars, 
primarily by delaying to higher densities, or even aborting the appearance of S _ hyperons, as 
shown in Fig. 4. 
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Figure 4. RMF calculations of baryon & lepton fractions in neutron star matter for different 
scenarios of hyperon nuclear potential depths Uy . Figure adapted from Ref. [19]. 



The G-matrices constructed from Nijmegen soft-core potential models have progressed 
throughout the years to produce X repulsion in symmetric nuclear matter, as demonstrated 
in Table 3 using the parametrization 

Vr = V »+\ V i T A-ts • (4) 



Table 3. Isoscalar and isovector £ nucleus potentials, Eq. (4) in MeV, calculated in Nijmegen 
soft-core potential models [20] at &p = 1.35 fm -1 , corresponding to nuclear-matter density. 





97f 


04d 


06d 


08a 


08b 


phenom. 


Ref. 




-13.9 


-26.0 


-1.2 


+13.4 


+20.3 


+30 ± 20 


[2, 16] 




-30.4 


+30.4 


+52.6 


+64.5 


+85.2 


» +80 


[21] 



20 



-20, 



- N 11™ al (k) with fss2 




k F [frrf'] k F [frrf'] 

1.45 1.30 

1.40 1.25 

1.35 1.20 



Uj al (k) with xEFT 




k F [fm 1 k F [fm 1 

— 1.15 0.9 

— 1.07 0.85 

— 1.0 0.75 



1 2 3 4 5 
k [frrf 1 ] 

(a) 



1 2 3 4 5 6 
k [fm -1 ] 



160 



120 

— 
E 

i 80 

5 

i 40 

Z 

< 



-40 



EN 1 S T=3/2 

fss2 

XEFT 

bare fss2 


"\ EN 3 S, T=3/2 
--__\bare fss2 

\\ 


''bare xEFT 


fss2 

XEFT 



2 3 
k [frrT 1 ] 





(b) 



2 3 
k [frrT 1 ] 



Figure 5. Left: Isoscalar £ nuclear potentials calculated in two models, Refs. [22, 23]. 
Right: (SA r ) T=3/ / 2 potentials calculated in these models. Figure adapted from Ref. [24]. 



In the latest Nijmegen ESC08 model [20], this repulsion is dominated by repulsion in the 
T = 3/2, 3 Si — 3 Z?i XiV channel where a strong short distance Pauli exclusion repulsion for 
quarks arises in SU(6) quark-model RGM [22] and in chiral EFT [23] calculations, as seen 
on the r.h.s. of Fig. 5. These model calculations also lead to £ nuclear repulsion, shown in 
momentum space on the l.h.s. of the figure. A strong repulsion appears also in a recent SU(3) 
chiral perturbation calculation [25] which yields Vq S ~ 60 MeV. Phenomenologically Vq S > 
and > 0, as listed in the table, and the resulting S-nuclear potential V§ is repulsive. 2 

4. K-nucleus potentials from K~ atoms 



The gross features of low-energy KN physics are encapsulated in the leading-order Tomozawa- 
Weinberg (TW) vector term of the chiral effective Lagrangian [6]. The Born approximation for 
the ^-nuclear potential due to this TW interaction term yields a sizable attraction: 
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for po = 0.16 fm -3 . Iterating the TW term plus the less significant NLO terms, within an 
in-medium coupled-channel approach constrained by the KN — ttT, — ttA data near the KN 
threshold, roughly doubles this ^T-nucleus attraction [29]. A major uncertainty in these chirally 
based studies arises from fitting the A(1405) resonance by the imaginary part of the (7rIl)/ = o 



2 In the case of |;He, the only known quasibound E hypernucleus [26, 27], the isovector term provides substantial 
attraction owing to the small value of A towards binding the T = 1/2 hypernuclear configuration, while the 
isoscalar repulsion reduces the quasibound level width [28]. 
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Figure 6. Left: Vr in Ni from global fits to 65 K~ atom data points, with x 2 values 
in parentheses [2]. The shaded area approximates the uncertainty in the Fourier-Bessel fit. 
Right: Overlap of K~ atomic 4/ radial wavefunctions R squared with Ni matter density p r , 



amplitude calculated within the same coupled channel chiral scheme. Yet, irrespective of this 
uncertainty, the A(1405) which may be viewed as a K~p quasibound state quickly dissolves in 
the nuclear medium at low density, so that the repulsive free-space scattering length a K - p , 
as function of p, becomes attractive well below pq. Adding the weakly density dependent 
1 = 1 attractive scattering length a^- n , the resulting in-medium KN isoscalar scattering length 
^o(p) = \{o-K-p{p) + a K-n{p)) translates into a strongly attractive V K [30, 31]: 

Vg(p) ~ -—Re bo(po) Po — « -HO - (MeV) . (6) 
Pkn Po Po 

Shallower potentials, Vg(po) ~ -(40 - 60) MeV, were obtained by imposing a Watson- 
like self-consistency requirement [30, 32]. It turns out, however, that stronger attraction, 

Vr(po) (80 - 90) MeV, arises in similar chiral approaches [33] when imposing the same 

requirement while considering the energy dependence of the in-medium KN scattering amplitude 
below threshold [34]. 

Comprehensive fits to the strong-interaction shifts and widths of AT~-atom levels, begun in 
the mid 1990s [35] , have yielded DD deeply attractive and strongly absorptive optical potentials 
with nuclear-matter depth -V^f(po) ~ (150 - 200) MeV at threshold [35]. The l.h.s. of Fig. 6 
illustrates for 58 Ni the real part of A'-nucleus potentials obtained from a global fit to the data in 
several models and, in parentheses, the corresponding values of x 2 for 65 AT~-atom data points. 
A model-independent Fourier-Bessel (FB) fit [36] is also shown, within an error band. Just 
three terms in the FB series, added to a tp potential, suffice to achieve a x 2 as l° w as 84 and to 
make the potential extremely deep, in agreement with the density-dependent best-fit potentials 
DD and F. In particular, potential F provides by far the best fit ever reported for any global 
AT~-atom data fit [37], and the lowest x 2 value as reached by the FB method. 

Shown on the r.h.s. of Fig. 6 are overlaps of the 4/ atomic radial wavefunction squared with 
the matter density p m in 58 Ni for two of the models exhibited on the l.h.s. of the figure. The 4/ 
atomic orbit is the last circular K~ atomic orbit from which the K~ meson undergoes nuclear 
absorption. The figure demonstrates that, whereas this overlap for the shallower tp potential 
peaks at nuclear density of order 10% of po, it peaks at about 60% of po for the deeper DD 
potential and has a secondary peak well inside the nucleus. The double-peak structure indicates 



Table 4. Full and reduced K atom data set fits. The reduced set consists of 2p, 3d, 4f, 5g, 7i 
shifts, widths and yields in C, Si, Ni, Sn and Pb targets, respectively [39]. 
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the existence of a K~ strong-interaction I = 3 quasibound state for the DD potential. It is 
clear that whereas within the tp potential there is no sensitivity to the interior of the nucleus, 
the opposite holds for the density dependent F potential which accesses regions of full nuclear 
density. This owes partly to the smaller imaginary part of F. 

Given the repercussions of deeply attractive potentials on the equation of state of dense 
matter, it is important to explore the stability of these best-fit solutions to variations in the 
data selection and the fitting procedure. The most obvious question to ask is whether the 
resulting best-fit potentials depend strongly on the size, composition and accuracy of the data 
set studied. Regarding size and composition, following an earlier discussion [38] it has been 
observed recently [39] that the DD deep potentials and the DI relatively shallow potentials, as 
well as the superiority of DD to DI in terms of quality of fit, persist upon decreasing the size of 
the data set. This is demonstrated in Table 4 upon reducing the 65 data point global set down 
to 15 data points from five targets spread over the entire periodic table (C, Si, Ni, Sn, Pb). 
Similar results hold for any four out of these five targets. It makes sense then to repeat some 
of the 30 — 40 years old K~ atom measurements, making use of modern techniques, in order to 
acquire a minimum size canonical set of data with reduced statistical errors and with common 
systematics. For the specific set proposed in Ref. [39], the (lower level) widths are directly 
measurable yet not excessively large to make it difficult to observe the feeding X-ray transition 
above the background. Similarly, the relative yields of the upper to lower level transitions are 
of the order of 10% and higher. Fitting to such a data set with improved accuracy could resolve 
the issue of deep vs. shallow potentials and determine how deep is 'deep'. 

A fairly new and independent evidence in favor of extremely deep ^T-nucleus potentials is 
provided by (K~ , n) and [K~ ,p) spectra taken at KEK on 12 C [40] and very recently also 
on 16 [41] at p K - = 1 GeV/c. The 12 C spectra are shown on the l.h.s. of Fig. 7, where 
the solid lines represent calculations (outlined in Ref. [43]) using potential depths in the range 
160 — 190 MeV. The dashed lines correspond to using relatively shallow potentials of depth 
about 60 MeV which may be considered excluded by these data. However, Magas et al. [44] 
have recently expressed concerns about protons of reactions other than those directly emanating 
in the (K ~ , p) reaction and which could explain part of the bound-state region of the measured 
spectrum without invoking a very deep ^-nuclear potential. A sufficiently deep potential would 
allow quasibound states bound by over 100 MeV, for which the major KN — > 7rE decay channel is 
blocked, resulting in relatively narrow K-nuclear states. Of course, a fairly sizable extrapolation 
is involved in this case using an energy-independent potential determined largely near threshold. 
Furthermore, the best-fit Vj K imaginary depths of 40 — 50 MeV imply that ^-nuclear quasibound 
states are broad, as studied in Refs. [37, 47]. 

A robust consequence of the sizable ^-nucleus attraction is that K~ condensation, when 
hyperon degrees of freedom are ignored, could occur in neutron star matter at about 3 times 
nuclear matter density, as shown on the r.h.s. of Fig. 7. Comparing it with Fig. 4 for neutron 
stars, but where strangeness materialized through hyperons, one may ask whether K mesons 
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Figure 7. Left: KEK-PS E548 missing mass (K~ ,n) (upper) & (K~~,p) (lower) spectra on 12 C 
at p K - = 1 GeV/c [40]. Right: calculated neutron-star population as a function of density. The 
neutron density stays nearly constant once kaons condense [42]. 




Figure 8. Saturation of K~ separation energy in RMF calculations of multi K~ nuclei as a 
function of the number k of K~ mesons [48]. Left: RMF model dependence. Right: dependence 
on RMF ingredients. 



condense also in the presence of hyperons. This question was posed within RMF calculations of 
neutron star matter long ago and answered negatively [45, 46], but only recently it was posed 



for strange hadronic matter in Ref. [48] by calculating multi-K nuclear configurations. Fig. 8 
demonstrates a remarkable saturation of K~ separation energies B K - calculated in mx\\t\-K~ 
nuclei, independently of the applied RMF model as shown on the l.h.s. for three different nuclear 
RMF schemes. The r.h.s. of the figure demonstrates that this saturation persists already in 
the most straightforward a + oj model, primarily owing to the repulsion induced by the vector 
oj field between like K mesons. The additional vector fields p and 4> only add repulsion, thus 
strengthening the saturation. The effect of Vj K is noticeable only below 100 MeV, as seen 

by the departure of the lowest green line with respect to the lowest red line. 

The saturation values of B K - do not allow conversion of hyperons to K mesons through the 
strong decays A — >■ p + K~ or H~ — >■ A + K~ in multi-strange hypernuclei, which therefore 
remain the lowest-energy configuration for multi-strange systems [49] . This provides a powerful 
argument against K condensation in the laboratory, under strong-interaction equilibrium 
conditions [48]. It does not apply to kaon condensation in neutron stars, where equilibrium 
configurations are determined by weak-interaction conditions. This work has been recently 
generalized to multi- K~ hypernuclei [50]. 



Appendix: density dependent optical potentials 



Here we specify the functional form of two density dependent optical potentials used in studies 
of hadronic atoms. For a recent application to K~ atoms, see Table 1 of Ref. [37]. 

• The DD form is based on modifying the effective scattering length bo (e.g. Eq. (6)): 

V*( r ) ~ b p(r) => b -> b + B {^} a , a > , (A.l) 

PhN P0 

where po = 0.16 fm~ 3 is a central nuclear density. It is possible then to respect the 'low 
density limit' by keeping bo fixed, bo = &o rcc , while varying the parameters Bq and a. 



• The F form is based on modifying bo as follows: 

b -> B F(r) + b [1 - F(r)} . (A.2) 
The density-like function F(r) is defined as 

F(r) = , x = . (A.3) 

e x + 1 a x 

Clearly, F(r) — > 1 for r « R x which defines an internal region and similarly [1 — F(r)] — > 1 
for r >> R x which defines an external region. Thus R x forms an approximate border 
between internal and external regions, and if R x is close to the nuclear surface then the two 
regions do correspond to the high-density and low-density regions of nuclei, respectively. 
In global fits across the periodic table, R x is parametrized as R x = RxoA 1 ^ + 5 X and the 
parameters Bq, R x o and S x are varied upon in the least-squares fit, while gridding on values 
of a x around 0.5 fm. The parameter bo may be held fixed at its free hN value, but the 
results often depend very little on its precise value. 
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